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Abstract 

The XX spin-chain with non-Hermitian diagonal boundary conditions is shown to 
be quasi-Hermitian for special values of the boundary parameters. This is proved by 
explicit construction of a new inner product employing a " quasi- fermion" algebra in mo- 
mentum space where creation and annihilation operators are not related via Hermitian 
conjugation. For a special example, when the boundary fields lie on the imaginary axis, 
we show the spectral equivalence of the quasi-Hermitian XX spin-chain with a non-local 
fermion model, where long range hopping of the particles occurs as the non-Hermitian 
boundary fields increase in strength. The corresponding Hamiltonian interpolates be- 
tween the open XX and the quantum group invariant XXZ model at the free fermion 
point. For an even number of sites the former is known to be related to a CFT with 
central charge c = 1, while the latter has been connected to a logarithmic CFT with 
central charge c = — 2. We discuss the underlying algebraic structures and show that 
for an odd number of sites the superalgebra symmetry [7(g[(l|l)) can be extended from 
the unit circle along the imaginary axis. We relate the vanishing of one of its central 
elements to the appearance of Jordan blocks in the Hamiltonian. 
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1 Introduction 



Recent years have seen increasing interest in non-Hermitian quantum Hamiltonians H and 
how to give them a physical sound interpretation. In this article we revisit an exactly- 
solvable, one-dimensional, discrete system, the XX-Hamiltonian with non-Hermitian diagonal 
boundary conditions, i.e. 

1 M_1 z 4- R z 

H =^H K<+i + + 2 — . «^6C. (1) 

Here a^ v,z denote the Pauli matrices acting in the m th copy of an Af-fold tensor product 
y®M Q £ a two-dimensional complex vector space V = Cv + © Cu_. For real values the 
parameters a, (3 € R have the physical interpretation of external magnetic fields located at 
the boundary. Here we are interested in investigating the case of complex fields, a, j3 £ C. 
Another variant of the XX-Hamiltonian we shall consider is 

1 M ~ 1 4- R M 

H '=2^ + < a rn+i ~ P < ~ « = - ^y 1 £ < ( 2 ) 

m— 1 m— 1 

which up to boundary terms coincides with the closed or periodic XX spin-chain Hamiltonian 
in an external magnetic field. Both Hamiltonians can be diagonalized via a Jordan- Wigner 
transformation for arbitrary complex values of the boundary parameters a, /?. However, the 
Hamiltonians H, H' with a, (3 € C are in general non-Hermitian, 

HjtH*=H(a,p), (3) 

and one needs to explain how a meaningful quantum mechanical system in terms of iJ, H' can 
be defined. On physical grounds the time evolution operator U(t) — exp(itH), t > ought 
to be unitary and, hence, H needs to be Hermitian. Thus, a non-Hermitian Hamiltonian, 
H H* , appears at first sight to be in contradiction with conventional quantum mechanics. 
However, under certain assumption the Hamiltonian might turn out to be quasi-Hermitian 
[TJ [5] • That is, there exists a positive, Hermitian and invertible operator n satisfying 

ijH = H*r/ . (4) 

This allows one to either introduce a new inner product, 

(v,w) n := {v,r}w), v,w € Sj (5) 

on the state space ft with respect to which H becomes Hermitian or perform a similarity 
transformation to a new, Hermitian Hamiltonian with respect to the original inner product, 

h = ^l 2 ^- 1 ' 2 . (6) 

Given a non-Hermitian Hamiltonian it is in general rather difficult to determine whether 
it is quasi-Hermitian. A special subclass of quasi-Hermitian systems where this turns out 
to be easier are those distinguished by PT-symmetry, i.e. Hamiltonian systems where the 
eigenvectors can be chosen such that they are eigenvectors under a joint parity and time 
reversal transformation; see e.g. [3] and references therein. In this case the quasi-Hermiticity 
operator 7/ enjoys further constraints which we will discuss in the text below. Employing the 
exact solvability of the non-Hermitian XX spin-chain we will establish for which values of 
the boundary parameters it is quasi-Hermitian and show its PT-symmetry. 
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1.1 Long range bulk interaction from non-Hermitian boundary fields 

The above non-Hermitian Hamiltonians are of special interest because of their underlying 
algebraic structure and we will show that it plays an important role in the interpretation of 
the new inner product. Setting a = —(3 = y— 1 the Hamilton ([1]) is the {7 9 (sZ2)-invariant 
XXZ spin-chain evaluated at the free fermion point q = y/—l [4j [5] . In the thermodynamic 
limit it has been suggested that the system is closely connected with critical dense polymers 
effectively described by a logarithmic conformal field theory with central charge c = — 2 
[7J 18]. Moreover, it has been stressed that such models should be described in terms of 
non-local degrees of freedom, such as "connectivities", see e.g. [Jj. 

On the other hand, setting a — (3 = the chain (fTJ) becomes Hermitian and is related to 
the Ashkin- Teller model, see [4] . In the thermodynamic limit the system is now described by 
an ordinary conformal field theory with central charge c = 1 . 

Thus, the non-Hermitian Hamiltonian 

^ M-l 

H 9 = t;J2 + «+i + - = H*_ g , < g < 1, (7) 

rn—l 

interpolates between these two special cases. We will show for small values of the coupling 
constant g its spectral equivalence with the Hermitian Hamiltonian 

M-n 

h g 2 = Px U) (9 2 ) [c* x c x+n - c x c* x+n ] , (8) 

n>0 x=l 

where c* , c x are fermionic creation and annihilation operators at lattice site x. The hopping 
probability between a site x and its n th neighbour is encoded in the real coefficients p x n ^ 
which only depend on g 2 and vanish for n even, j4 2 ™' =0. At g = we only have nearest 
neighbour hopping, 

AI-l 

ho = - ( c *x c x+i - c x c* x+1 ) . (9) 

x = l 

Our perturbative calculation will show that as g increases the bulk interaction becomes more 
and more long-range by successively "switching on" the various coefficients p x 2n+1 ^ starting 
from the boundary sites x — 1, M . More precisely, we find up to order g s that the nonvan- 
ishing contributions are 

(1) 128. 9 2 + 8ff 4 + ff 6 , t 

Px = 1 (Sx.l+Sx.M-l) 

_ Sg+3g_^ s + s _^ + L(j + $ xM _ 3 ) + o(g s ) 
512 25o 

P * ] = 20g4 25 + 6 3g6 + **,M-3) + ^(Sx,2 + 6 X , M -,) + 0(g s ) 

P X 5) = -^(6x,i+5 x ,M- 5 )+0(g s ) (10) 
512 

Figure 1 indicates which hopping amplitudes, i.e. the probabilities for a fermion to jump 
between two lattice sites, arc modified up to order g 8 . We have only depicted the additional 
contributions for g > 0: whenever an arc connects two lattice sites x and x + n then there is 
a non- vanishing coefficient p x n ^ . 



2 



12 3 4 5 6 



M-5M-4 M-3 M-2 M-l M. 



Figure 1. Graphical depiction of the long range hopping. 

Both variants ([7]) and §8§ of the Hamiltonian have their advantages: 

• The Hermitian Hamiltonian h serves physical intuition. It can be directly interpreted 
and clearly shows the long range nature of the interaction which is not apparent in ([7]). 
It suggests that the long range nature of the interaction with increasing g is behind the 
singular change from an "ordinary" conformal field theory with c = 1 for < g < 1 to 
a logarithmic one with c = — 2 at g = 1. 

• In contrast, the variant ([7]) highlights the algebraic properties, such as integrability and 
quantum group invariance. The PT-symmetry and quasi-Hermiticity of the quantum 
group invariant XX Z spin-chain for higher roots of unity q = exp(i7r/r), r > 2 have 
been discussed in another paper [9]. 

Note that the appearance of a long range bulk interaction does not contradict physical 
intuition. Naively one would argue that boundary terms will become unimportant in the 
thermodynamic limit when the number of sites tends to infinity. However, the applicability 
of this statement crucially depends on the nature of the boundary conditions, see for instance 
[lOj . In the present article we consider a different, novel example of this phenomenon. The 
complex boundary terms render the Hamiltonian non-Hermitian, hence they imply a drastic 
non-local change in the mathematical structure of the state space: the introduction of a new 
inner product. We will see that this is the case by observing that the similarity transformation 
which maps H g into h g 2 is highly non-local. 

1.2 Outline of the article and summary of results 

The purpose of the present article is to relate the notions of PT-symmetry and quasi- 
Hermiticity to another simple but non-trivial model. While the XX spin-chain with non- 
Hermitian boundary terms has been discussed previously in the literature (see e.g. jTTJ Q21 
[T3l [14j [15] ) it has not been investigated for which parameter values a,/3 the Hamiltonian is 
quasi-Hermitian and what the corresponding Hermitian systems are. We will show that for 
a = j3 inside the unit disc the Hamiltonians ([I]), (|2|) are quasi-Hermitian, this condition on 
a, (3 ensures also PT-symmetry. We then will discuss for which values the mentioned Hamil- 
tonians cease to be quasi-Hermitian and when their spectrum contains complex eigenvalues. 
The special case a = (3 of the non-Hermitian XX spin-chain, which includes ([7]), has not been 
discussed in detail previously. 

For the benefit of the reader we summarize the new results contained in this article: 

• After a brief review of the notions of quasi-Hermiticity and PT-symmetry on the lattice 
in Section 2, we establish in Section 3 that the spectra of the Hamiltonians |T|) and 
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@ are real for boundary parameters within the unit disc, a = (3 £ C and \a\ < 1. 
Performing a Jordan- Wigner transformation the XX spin-chain can be reformulated 
as a non-trivial fermion model and its spectrum can be described in terms of quasi- 
particle excitations in momentum space. We will give an elementary proof that all quasi- 
momenta (Bethe roots) lie on the unit circle and state for a — —0 — ig, < g < 1 a set 
of palindromic polynomials whose roots give the quasi-momenta k and the associated 
energies 2 cos k. 

• In Section 4 we highlight that unlike in the case of real boundary fields a, (3 £ M, 
the Jordan- Wigner transformation does not lead to a well-defined fermion algebra in 
momentum space, but instead one has two sets {c* k ,Ck\ and {d* k ,dk} of creation and 
annihilation operators which satisfy the relations 

[c*,4']+ = 8k,k', [h, 4']+ = = 0, d k ■ (11) 

Here k = — i In z £ R is a quasi- momentum, [A, B] + = AB + BA is the anti-commutator 
and * denotes the Hermitian adjoint with respect to the original inner product, where 
the Hamiltonian is non-Hermitian. Employing the above operators we explicitly con- 
struct an i] which not only renders the Hamiltonian Hermitian but obeys the more 
restrictive condition 

ric* k = d* h r) and i]d k = c k r] . (12) 

That is, with respect to the new inner product ([5]) {c* k ,dk\ satisfy the canonical anti- 
commutation relations and are the Hermitian adjoint of each other. 

• Also in Section 4 we extend for a — —(3 — ig, < g < 1 the well known U q —i{sl2)- 
symmetry of ([7]) from the unit circle, g = 1, into the unit disc, < g < 1 and 
discuss the invariance of the Hamiltonian for odd and even number of sites. Both 
cases show remarkable differences. The quantum group symmetry is closely connected 
with so-called fermionic zero modes and we show that at the coupling values g where 
the anticommutator of the associated fermionic creation and annihilation operators 
vanishes the Hamiltonian possesses non-trivial Jordan blocks. Thus, the Hamiltonian 
ceases to be quasi-Hcrmitian. For even numbers of sites this happens at g = 1, i.e. on 
the unit circle, as was observed previously [TJ. Here we show that it also happens for 
odd numbers of sites albeit at the value g max = y (M + 1)/(M — 1) which approaches 
the value one in the limit M — » oo. It also signals the onset of complex eigenvalues for 
values of g > (? ma x, i-e. outside the unit disc. 

• In Section 5 we present a perturbative calculation of rj in terms of the coupling pa- 
rameter < g <C 1 in 0. While we follow closely the steps previously put forward in 
the literature [161 117) , we give a novel derivation of the coefficients in the perturbation 
series expansion. Using these results we obtain the expression {SJ. We also will derive 
some closed expressions for small number of sites M = 3, 4 and 5 in Section 6. For g = 1 
and odd lattice sites we formulate a conjecture on an alternative way of computing the 
new inner product using the Temperley-Lieb algebra. 

• In Section 7 we will investigate briefly the case of general complex boundary fields 
inside the unit disc, i.e. when a = f) = ge %e with 8 ^ it/ 2 and real in order to make 
contact with the discussion in [IT]. 

• Section 8 contains the conclusions. 

2 Quasi-Hermiticity and PT-invariance 

The state space of the non-Hermitian systems (JJ and ^ is a spin-chain of M sites repre- 
sented in terms of the tensor product V® M of the two-dimensional complex vector space V 
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with orthonormal basis {v±} such that a z v± = ±v±. The tensor product is then spanned 
by the vectors 

{|ei,...,ejvf) = v Bl ®---®«e M : £m = ±l} (13) 

This particular choice of basis vectors is motivated by the axial symmetry of the Hamiltonians 
H, H' which commute with the total spin operator, 

1 M 

S Z = ^E<> [H,S*} = [H',S z ]=0. (14) 

z m=l 

The matrices (fTJ) and ([2]) are defined with respect to the inner product 

M 

(ei,...,£ M Wi,-,£'m) = Y[( v ei,V E >), (v e ,V e i) = S ee ' ■ (15) 

i=l 

As already pointed out in the introduction both Hamiltonians, H and H 7 , are in general 
non-Hermitian for arbitrary complex boundary parameters a, 0. However, in the above basis 
(|13[) both of them are symmetric, 

H = H t and H' = H n . 

We now wish to determine for which values there exists a map r\ : V® M — > V® M possessing 
the following properties: 

1. rj is Hermitian, 77 = rf , invertiblc, dot 77 > 0, and positive definite, r/ > 0. 

2. rj intertwines the Hamiltonian H (respectively H') with its Hermitian adjoint, 

77^ = ^*77. (16) 

Because of the axial symmetry present in H and H' we add the following requirement 
which implies that 77 intertwines H and H' at the same time, 

[r),S z ]=0. (17) 

Provided that such a map r\ exists we define a new inner product and with it a new Hilbert 
space structure via 

(v,w) v :=(v,r)w), v,w eV® M . (18) 

Using the intertwining property (2) it is obvious that the Hamiltonian is Hermitian with 
respect to the new inner product. The properties listed under (1) ensure that the new inner 
product is well-defined. 

There is an alternative to introducing a new Hilbert space structure. Since rj > it 
follows that there exists a unique positive definite square root 772 > which allows one to 
define the Hermitian Hamiltonian 

h = rj^Hr]~^ (19) 

with respect to the original inner product. 

Both approaches have their advantages and disadvantages. At first glance one might 
prefer to work in the original Hilbert space with the "gauge transformed" Hamiltonian h, 
as this allows for a direct physical interpretation without having to insert another operator 
when taking scalar products. However, for practical purposes it is often more feasible to work 
with H and 77. Even if it is possible to construct 77 explicitly, the computation of its square 
root and, thus, the calculation of h is another technically complicated step. We will perform 
this computation for ([7]) to low orders in g employing perturbation theory in Section 5. In 
contrast, the original Hamiltonian H has usually a simpler form which allows one to read off 
certain symmetries. The XX-chain with imaginary boundary fields is a concrete example. 
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2.1 Parity, Time and Spin Reversal on the lattice 

On the set of basis vectors (fTBf we define the parity reversal operator P by 

P\e 1 ,...,e M ) = \e M ,£M-i,—,si) (20) 

and extend its action to the whole space by linearity. From this definition it is immediate to 
see that 

PHP = H{f3,a). (21) 
The time reversal operator T acts on the basis vectors as identity, 

T|ei,..,£ M ) =|£i,...,£m>, (22) 

but is defined to be antilinear, whence any matrix (such as the Hamiltonian) is transformed 
into its complex conjugate under the adjoint action of T, 

THT = H = H(a,f3) = H* . (23) 

Below we will define discrete wave functions which will be transformed into their complex 
conjugates under the operator T . This justifies the identification of T with time-reversal. 
Upon imposing the constraint 

a = (24) 
the Hamiltonian turns out to be PT-invariant, 

a = $ : [PT, H] = . (25) 

Note, however, that PT-invariance of the Hamiltonian is not a sufficient criterion to 
ensure real eigenvalues, this only follows once it is established that the eigenvectors of the 
Hamiltonian can be chosen to be simultaneously eigenvectors of the PT-operator. This is not 
automatically implied by PT-invariance of the Hamiltonian as time reversal T is an antilinear 
operator; see e.g. [18]. 

For later purposes we also discuss the behaviour under the spin-reversal operator, 

M 

RH(a,(3)R = H(-a,-(3), R= \[ <j* m . (26) 

m—l 

Thus, in general spin-reversal symmetry is broken in the presence of boundary fields. 

Provided that the Hamiltonian is not only PT-invariant but satisfies the slightly more 
restrictive constraint 

PHP = THT = H*, (27) 

as it is the case here for a = /3, it is natural to impose further constraints on rj. Namely, we 
wish to have 

P V P = Tr/T = if 1 . (28) 

These conditions are compatible with the intertwining property, e.g. 

Pj]HP = Pr]P H* = H P?/P = PH*r]P . 

In essence we are demanding that parity-reversal gives the same rj up to inversion. The same 
applies for time-reversal. 

Often it is beneficial to define an additional operator C introduced by Bender and col- 
laborators (for references see [3]) by setting 

C := Pry . (29) 

This operator might turn out to have a simpler expression than rj itself. The aforementioned 
properties of rj imply the identities 

C 2 = 1, [PT, C] = 0, [H, C] = . (30) 

The C-operator turned out to have an elegant algebraic expression for the quantum group 
invariant XXZ spin-chain; see [9]. 
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3 Spectrum and eigenvectors of the Hamiltonian 

As mentioned in the introduction the Hamiltonians H, H' can be diagonalised in terms of 
free fcrmions. Using the well known Jordan- Wigner identities 

c * = ( II o* v J c* = ( n a* v J <rt, and n x = c\c x = (31) 

\y<x J \y<x ) z 

we introduce fermion creation and annihilation operators in "position space" satisfying the 
canonical anti-commutation relations (CAR), 

[c x ,c y } + = [c*,c*] + = and [c*, Cy]+ = 5 XjV , x, y = 1, 2, M . (32) 

For later purposes we also state the transformation properties under parity, time and spin- 
reversal, 

Pc x = c M +i- x P, Tc x = c x T, Rc x = (-l) x+1 c x R. (33) 
In terms of the fermion algebra the Hamiltonians can be rewritten as 

M-i a + (3 
H = - ^ [ c *x c x+i ~ c xC* x+1 ] + a m + (3 n M ^— , (34) 



and 

M - 1 r i a 
H = — y c x c x+ i - c x c x+1 + p n x + a n x+1 — 



(35) 



Below we shall constrain the boundary parameters a, (3. For the moment we leave them ar- 
bitrary. We now introduce a "discrete wave function" i\) z depending on a complex parameter 
z G C by defining 

M 

K = 22ip z (x;a,(3)c*, (36) 

x=l 

with 

ip z (x;a,(3) = z x - A(z;a,0) z~ x . (37) 

This ansatz is physically motivated: it is a superposition of two (discrete) plane waves, one 
incoming and one reflected. Here z — exp(ik) with k being the quasi-momentum and A is 
the reflection coefficient. Note that we allow here for complex momenta k. 
Employing the canonical anticommutation relations one easily finds that 

c x H = Hc* x + c* x _ x + c* +1 , l<x<M, 
c[H = H c j — ac\ + c* 2 , 
c mH = H°M + c M-i ~ P C M ■ 

Therefore, one has 

[H,c z ] = -(z + z' 1 )c* z (38) 
provided the coefficient A in the wave function ([3T)l obeys the identities [1] 

1 + a/z /3 + z- 1 y ' 

As both equations have to hold simultaneously this imposes a constraint on the allowed 
values for the parameter z which are specified as the roots of a polynomial equation of 
order 2M + 2. From (13"5)) we infer that the spectrum of the Hamiltonian is composed of 
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quasi-particle excitations with energy e — z + z 1 by successively acting with (1361) on the 
pseudo- vacuum vector |0) = V-i ® • • • ® U-i, 

# zi) = (-^A E (* + zt 1 )) ■••> *«> > (4°) 



where 

...,«,) =^...^,|0) , . (41) 

Thus, it depends on the nature of the solutions z of the equation (|39|) whether the spectrum 
of the Hamiltonian is real. 

Proposition 1 Let a = in {7p and assume that a is inside the closed unit disc, i.e. \a\ < 1. 
Then the solutions of H39\) all lie on the unit circle, hence the quasi-momenta kj = —i In Zj 
are real. 

Proof. Let us rewrite the Bethe ansatz equations (|39|) as 



zM _z + a_ = z _ M 

za+l z- l a + l 

and define the maps f(z) — z M (z + a)/(za + 1) and g(z) = f(z~ 1 ). Notice that the above 
Bethe ansatz equations are invariant under complex conjugation and z — > z -1 , which reflects 
the PT-invariance of the Hamiltonian. Thus, we can assume without loss of generality that 
there exists a solution z with \z Q \ < 1. For \a\ < 1 the image of the closed unit disc under 
the map / lies again in the closed unit disc. In contrast the image of the closed unit disc 
under the map g lies outside of the open disc. Hence, any solution to (|4"2)l with \a\ < 1 must 
lie on the boundary of these two image regions, i.e. the unit circle. H 

Remark. Motivated by the previous proposition we shall henceforth use the parametri- 
sation 

a = = ge ie , g > (43) 

for the boundary parameters. The Hamiltonian therefore depends now only on two real 
parameters, g > and < 9 < 2tt. 

Note that there are only M relevant solutions to the Bethe ansatz equations. The trivial 
roots z = ±1 do not occur in the spectrum of the Hamiltonian which allows one to reduce 
the problem of solving (|39p to finding the roots of the following palindromic or self-reciprocal 
polynomial 

M-l A/-1 
f(z) = z 2M f(z- x ) = z 2M + 1 + (1 + g 2 ) z2m + 2 9 cos0 z2m+1 ■ ( 44 ) 

m—l rn—0 

Thus, all roots Zi occur in reciprocal pairs and there exists a unique polynomial 

M 

F{e)^\{{e~e l ) (45) 

i=l 

of order M whose roots are given by the single-particle energies (compare with fl40[) ) 

e l = z l + zT 1 . (46) 

Note that the ambiguity in the definition of the Bethe root Zj does not matter, as the 
wavefunction (|37| is simply rescaled by changing Zi — > z±. Thus, the problem of computing 
the eigenvectors and spectrum of the Hamiltonian is reduced to finding the roots £j of a 
polynomial F of degree M. For example, setting M — 8 and p = a + (3, a = 1 + a/3 we 
obtain 

F(e) = e s - e 6 (8 - a) + 5e 4 (4 - a) - 2e 2 (8 - 3a) + 2 - a - Ape + 10pe 3 - 6pe 5 + pe 7 . (47) 

We now specialize to the case of particular interest, a = —(3 = ig, j€ R, and present general 
expressions for the reduced polynomial F for all M . 
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3.1 Palindromic polynomials for purely imaginary boundary fields 

For the special choice 6 = ±7r/2 the palindromic polynomial (|4"4")) simplifies since all terms 
involving odd powers disappear. One easily convinces oneself, that this leads to the further 
simplification that all roots occur in pairs ±£j. Furthermore, if the number of sites M is odd 
one easily verifies that one has the roots Z{ = ±y/—l. The latter give rise to a "zero mode", 
£i =0, of the Hamiltonian 

where the corresponding wave function is given by 



sin ^ — ig cos ^ 



M+l - M=±r,2 

2 2 y 



A similar expression holds for z = — i. As we will see below this solution for the discrete wave 
function is connected with a [/(g[(l|l))-symmetry of the Hamiltonian. Dividing out the zero 
mode, we end up with a palindromic polynomial of degree 2(M — 1), 



Setting therefore 



we write once more 



4k+2 



M- 1 



fc=0 



k=0 



M-l 



, . M odd 

^, M even 



f(z) = Y[(z 2 -EiZ + 1), e % 



(48) 



(49) 



for the reduced polynomial. From the two alternative expressions for / we obtain a linear 
system of equations for the elementary symmetric polynomials e& = ek(£i, £ m ) in the roots 
Si. Solving this system we then define the polynomial 



2 m 



F(e) = ^(-l) fe efe e m - fe = H(e e t )(s + £i ) 



k=0 



i=l 



We now explicitly state this polynomial F for a = —f3 on the imaginary axis. We have to 
distinguish the cases of odd and even sites: 

M = 2m + 1. 



F(e) 



-) m E 



k=0 



{-) k e 2k 
(2fc+l)! 



(m + l + fc)! 2 (m + fc)! 



(m-fc)! 



(m — 1 — k)\ 



= 2 



sin [(M + 1) arccos |] + g 2 sin \(M — 1) arccos |] 
e\/4 - e 2 



(50) 
(51) 



M — 2m. 



F(s) 



(-) m E 



kJlk 



k=0 



{-Ye 
(2k)\ 



(m + k)l 2 (m-l + fc)! 
(to — fc)! ^ (m— 1 — fc)! 



= 2 



sin [(M + 1) arccos |] + g 2 sin [(M — 1) arccos |] 



(52) 
(53) 
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The expressions in terms of the inverse function arccos can be checked by rewriting the 
Bethe ansatz equations (j3T))) . At a = —j3 = ig the latter simplify to the transcendental 
equation, 

_ 2M (3 + z 1 + az- 1 2 m z 2 + g 2 2 _ sin[(M + l)C] _ 

+ z- 1 1 + az ~ 1 + zV sin[(M-l)C]' ~ Pl U ' 

(54) 

This is in agreement with the above expressions for the single particle energies as zeroes of 
(fSTj) and (|53|) . Note however that the polynomial expressions (|50|) . (|52|) are of advantage 
in the numerical computation of the single particle energies. From (|50|) one also easily 
spots the occurrence of another zero mode e — at g 2 = (M + 1)/(M — 1) for M odd 
which we will connect below with a representation of the universal enveloping algebra of the 
Lie superalgebra g[(l|l). We will also exploit the above polynomials and their expressions 
in terms of inverse trigonometric functions to derive approximations for the single particle 
energies in the interval < g < 1. 

3.2 Groundstate eigenvalues and central charges 

Using the results from the previous section we can expand the single particle energies as 
power series in the coupling parameter g 2 . Employing the equation 

F(e) = 

we find the following approximated expressions in the vicinity of the points g 2 = and g 2 = 1 
where the exact solutions are known: 

Trfc „ sin jf^r sin ■ 271 L 



Q<g<l: e k ^2co Sj ^-2g 2 M±±_M±L + (g*), k = l,...,M (55) 

and 

7rfc , 9N sinfytan^l „ „ r s 

0« 5 <l: £fc = 2cos — + (l-.g 2 ) M M +0((l-g 2 ) 2 ), fc = 1, ...,M— 1 . (56) 

If M is even we have to omit the value fc = M/2 where the above approximation is not valid 
since we then have a zero mode at g = 1 and the equation F(0) = becomes trivial. For M 
even we thus approximate the two missing eigenvalues which converge to e — at g = 1 by 



e± , = ±2 



2(1 - g 



2 \ 



M /2 y M(M + 2 — g 2 (M — 2)) ' 

This approximation is found by expanding the Bethe ansatz equations (153j) at e = 0. An 
example, M = 8, for these approximations is shown in Figures 2 and 3. The solid lines 
indicate the exact single particle energies S{ and the dashed lines the approximations. In the 



10 



vicinity of g = the approximations are depicted in Figure 2. 

approximated single particle energies e, for M = 8 



1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 




0.0 



0.2 



0.4 0.6 

S 



0.8 



1.0 



Figure 2. Approximations of the single particle energies in the vicinity of g = 0. 



For the approximations in the vicinity of g — 1 see Figure 3. 

approximated single particle energies e, for M = 8 
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-1.5 
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Figure 3. Approximations of the single particle energies in the vicinity of g = 1. 

Having found approximated expressions for the single particle energies £j, we can derive 
approximations for the groundstate eigenvalue Eq of ([7]) and discuss the finite size corrections 
of the system as it approaches the two points 5 = and g = 1. For the latter values the 
result is well known, see Table 1. 
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.9 = 


.9 = 1 


M odd 


p _ -, 2(M+1) w n / 1 v 


£J = 1 - 2^ - ^ + O(j^) 


M even 




E a = 1 - + + 0(^j) 



Table 1. Expansion of the groundstate energy with respect to the lattice size M. 
Using these expansions in the system size M one obtains via the general formula [33l [34] 

So = 2M/ 00 +/ s -^+0(4) (57) 



the effective central charge c e // = c — 12d min of the conformal field theory describing the 
system in the thermodynamic limit. Here d m i n is the smallest scaling dimension occurring in 
the theory and f^, f s are the bulk and surface free energy, respectively. One recovers the 
familiar central charges at g = 1 mentioned in the introduction: for an odd number of sites 
the conformal anomaly is c = 1 (d m - ln — 0) while it is c = —2 for M even (cZ m ; n = 0). Note 
that for M odd at g = we also formally obtain c e ff = —2, it appears however that this 
case has not been investigated further in the literature [35] . 

Let us now consider the approximate expressions for < g < 1. Summing the above 
expressions for the single particle energies we arrive at 



M odd, < g < 1 : 

Eq = — 
= 1 

M even, < g < 1 : 



M-l 

2 

£ 

fc=i 



1 



2(M+1) 



2 cos 



M+l 



sin • 



2(M+1) 
2(M+1) 7T 

tt + 6M 



(M + 1) sin 



2(M+1) ^ 2 



377 



2(M+1) 

+ 0(M- 2 ) + 0( 3 4 ) 



(58) 



E = - £ k = 1 - 



1 



1 



k=l 

4g 2 
+ 37 



2(M+1) 

2(M + 1) 7 



cos 



M+l 



(M - 



1 ) S111 2(i/+l) 



0(5 4 ) 



12M 



+ 0(Af- 2 ) + 0(g 4 ) 



(59) 



Thus in both cases we only see a change in the surface term f s up to order g 2 in the vicinity 
of .g = 0. A computation of the second derivative e'l(g 2 = 0) using the equation F(e) = 
shows that also the g 4 contribution to the groundstate only contributes to the surface energy. 
The computation of higher orders becomes cumbersome and it might be preferrable to rely 
on field theoretic methods instead, we briefly comment on this in the conclusions. 

Let us now turn to the case when g is the vicinity of g — 1 to see whether we encounter 
a shift in the central charge here. Again we sum the approximate expressions of the single 
particle energies to obtain the expression for the groundstate eigenvalue. Unlike for g = 
we cannot find exact expressions at this point. Instead we split the sum over the correction 
terms for g < 1 into two parts, 



^ sin % tan # 



M 



M 



k 



M 



1 \ - 1 1 x - 7rfc 
\ — \ cos — 

M 4^ cos # M 4" M 

k Ivl A; 



and for M > 4 employ the Euler-Maclaurin formula to obtain the following approximations 
for the first sum, 



M odd, < 



1 

< 1 : — 



•31n± 



1 



M ^ cos^ 



k=l 



M 



3tt 



— h 0( ) 

2M y M 2 > 



(60) 



12 



and 



Mi 

1 ^ 1 7- 12 In 4 +12 In 2 i , i 

Meven, 0< g < 1 : — V r r> M +0 — . 61 

M f- 1 cos w 12tt 2Af y M 2 ' V 7 

For M even we have another contribution from the energy level which becomes a zero mode 
at .9 = 1. Keeping <5 = 1 — g 2 fixed and choosing M large enough such that SM > 1 we find 
the asymptotic expansion 



2V2 2\/2(2-<5) 1 



^ - ir — jm^ + 5 := 1 - - 9 (62) 

However, for our purposes this is not a good approximation since it should vanish at g = 1. 
If we therefore expand first with respect to S we find instead, 



[25 (M-2) ,3/2 , ^,5/2, 

EM/2 = V 17- 7 /77TT ^ +°(' 5 ) 



(63) 



From this approximation we can not infer the correct finite size scaling properties. In order 
to extract conclusive results about the finite size scaling behaviour we therefore need to find 
the exact solutions for the single particle energies first. For now we have to leave this problem 
open, but we hope to address it in future work by field theoretic methods and exploiting the 
underlying algebraic structures which we highlight next. 



4 Quasi-Fermions and £7(0t(l|l))-invariance 

As we saw in the previous section the spectrum of the Hamiltonian is composed out of 
single particle excitations with energies £j = + z~ x which are "created" from a pseudo- 
vacuum applying the creation operator (|36j) . Naturally, we wish to find the corresponding 
annihilation operator in order to set up a Fermion or CAR algebra in quasi-momentum space 
which diagonalises the Hamiltonian. As the Hamiltonian is non-Hermitian the annihilation 
operator of momentum k — — ihxz cannot simply be the Hermitian adjoint of (|36p with 
respect to the original inner product. Instead we introduce the time-reversed annihilation 
operator 

M 

d z := Tc z T = } j ip z (x; a,0 = a)c x , (64) 

x=l 

where the conjugation with the time reversal operator is motivated by observing that 

THT = H * for (3 = a . 
Obviously, we then obtain the desired commutation relation 

[H,d s ] = (z + z- 1 )d M , (65) 

i.e. d z annihilates a quasi-particle of energy e = z + z^ 1 . 

Alternatively, we could also have employed the parity operator P in the definition (|64p . 
since PHP = H* . However, up to a phase factor this would lead to the same result due to 
the following identities of the wave function ij) z (x; a) := ip z (x; a, (3 = a), 

ip z (M + l-x;a) = -z M+1 ^ z (z; a) = -z M+1 ip z (x; a) . (66) 

Here we have assumed that z is a solution of (|39p and hence, lies on the unit circle. From 
(f6"6")) together with (f3"3"]) we infer that 

PTc* z PT = -z~ M ~ l c* z and Pc z P = -z~ M ~ l d z . (67) 



13 



This in particular implies that the system possesses PT-symmetry, i.e. the eigenstates of the 
Hamiltonian can be chosen to be eigenvectors of the PT-operator. However, this fact does 
not guarantee quasi-Hermiticity which we show next. 

Our strategy for establishing quasi-Hermiticity is to turn the operators (|36|) and (|64p 
into a well defined representation of the fermionic oscillator or CAR algebra. Under an 
appropriate renormalization of the wavefunction (|37|. 



V> z (x;a)-> a z ^ z (x;a) , = 1 

one verifies that the following anticommutation relations hold, 

[K x ,c*z 2 }+ = [4 1? 4 2 ]+ = and [c* Zi ,d Z2 } + = S ZUZ2 . (69) 



Here z\ , Zi are two solutions to the Bethe ansatz equations (j39[) . The anticommutation 
relations (|69|) are a direct consequence of the wave function identities 

^i } Zk {x)i) Zl {x)=8 k i and ^ i>z h (x)ip Zk (y) = S x , y , (70) 

x k 

where the index k labels the solutions to the Bethe ansatz equations. Notice that these 
identities lead to the inversion formulae 



k k 

and 

Cx=Y^ ^z k (x)c Zk = ^z k 0)4* ■ (72) 

k k 

While we have obtained the correct anticommuation relations in quasi-momentum space, the 
algebras generated by {c*,d z } do not possess the right ^-structure (anti- involution), where 
creation and annihilation operators are related by Hermitian conjugation. However, we now 
introduce the quasi-Hermiticity operator and the associated inner product with respect to 
which ([36| and (|64|) possess the right ^-structure. 

Theorem 2 The Hamiltonians {IJ), 0) are quasi- Hermitian for a — [3 and \ot\ < 1. Let 
Zj — expikj, j — 1,...,M be the M roots of the polynomial [44]) n( >t containing a reciprocal 
pair (kj ^ —kj> for all j,j'), then 

, - M 

H= I_-^2cosfc J :c* k .d kj . (73) 

3=1 

A similar expression holds for 0). The quasi-Hermiticity operator r\ : y® M — > y® M an d n s 
inverse are given by 

M 

»? = £ £ 4 3 . 1 ---4jo><o|4,„---4 1 (74) 

n=0 <— <»in 

and 

M 

T 1 = E E • • • fi L 1°) (°l%n ' ' ' ■ ( 75 ) 

n=0 <"-<™J« 

Besides quasi-Hermiticity of the Hamiltonian the map r\ also induces the conventional *- 
structure of free fermions in quasi-momentum space via the identities 

f]c\ = d* k rj and r/dk = Cfc?7, k = 1, M . (76) 
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Proof. First we notice that the creation operators (|36|) evaluated at the solutions z — Zi 
of the equations |39|) provide us with a basis in the state space, see (|4Tj) . Due to the 
symmetry of the equations (|39[) under the replacement z — > z , there are only M relevant 
distinct solutions, despite the fact that the polynomial order is 2M + 2. According to the 
anticommutation relations ([69)) these solutions then yield 2 M = dim V® M eigenvectors. The 
assertions then follow from the previous proposition showing that all Bethe roots lie on the 
unit circle for a = j3 and \a\ < 1 as well as employing (|40|) and (|69|) . H 

Corollary 3 The C-operator C — Prj has the following expression in terms of creation and 
annihilation operators, 

M 

C=Y, £ ( _ ) n eI (M + l)fe 1+ ..4A.J^....,J 0)(0| ^...^ (?7) 

n— 31 < -'-- < ^3n 

Thus, with respect to the quasi-particle basis kj n ) = c* k . |0)}^L ^ e C- 

operator is simply given in terms of the total quasi-momentum operator P = kj r and 
the quasi-particle number N = y\ c k _ d kj as 

C = (-1) . (78) 

Proof. An immediate consequence from (|oT|) . ■ 

Let us summarize the result: the fermion or CAR algebra {c x ,c*}^f—i in position space 
with respect to the original inner product is replaced by a CAR algebra {dk, c* k } k 4 =1 in quasi- 
momentum space with respect to the 77-product, 

(c* k v,w). n = (v,c k i]w) = {v,d k w)r, (79) 

and 

{dkV,w) v = (v,d* k -qw) = (v,c* k w) 71 , v,w G V® M . (80) 

Note that there is another copy of this CAR algebra, namely {d k , Cfc}^f =1 but with respect to 
the ?7" 1 -product. 

We have excluded the unit circle from the allowed range of the boundary parameter a. 
While the spectrum of H is also real in this case, it does not necessarily imply the existence 
of a positive definite operator 77, i.e. quasi-Hcrmiticity is a stronger condition. This is linked 
to the appearance of Jordan blocks in the Jordan normal form of the Hamiltonian as we will 
discuss now for the special case when a — ~(3 — ig, g > 0. Without loss of generality we 
can restrict ourselves to positive values of g, since g and —g are related via spin-reversal. We 
will show that the appearance of non-trivial Jordan blocks coincides with the vanishing of a 
central element in a representation of a specific superalgebra. 

4.1 "Deformed" quantum group symmetry on the imaginary axis 

Since the Hamiltonians |T]), |2]) and {7} can be expressed in terms of free fermions as we have 
just seen, there are several symmetries and algebraic structures associated with them. These 
have been previously investigated for a, (3 on the unit circle. Here we shall show that these 
symmetries and algebras can be extended inside the unit disc along the imaginary axis setting 
a = —[3 — ig, g > 0. We shall also relate the new inner product to the representations of 
these algebras. We start by reviewing the known symmetries and algebras for the case g = X, 
i.e. a, (3 on the unit circle. 

It is well-known that for g = 1 the Hamiltonian (|7|) possesses a J7 g (s^2)-symmetry with 
q = i = y/— 1 [5 . Below we shall see that for an odd number of sites, M 6 2N + 1, this 
symmetry can be extended along the imaginary axis for < g < \J (M + l)/(Af — 1). We 
also comment on a deformation of the Temperley-Lieb algebra [HI [55] (see [53] for a text 
book), which is the commutant of U q (sl2). First we recall the basic algebraic definitions. 
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Definition 1 The q-deformed enveloping algebra U q (sl2) is defined in terms of the Chevalley 
generators {E, F, K ±x } and the relations 

KR- 1 = K~ X K = 1, KE = q 2 EK, KF = q~ 2 FK, [E, F] = ^— ^ . (81) 

Definition 2 The Temperley-Lieb algebra TLm(<?) is the associative, unital algebra gener- 
ted by {e^fL^ 1 subject to the identities 



in 



„2 / , -1 



(q + q )ei, eiei±iei = ej, ae 3 = ejBi, \i - j\ > 1 . (82) 

Setting the deformation parameter to q = \A-T we obtain at g = 1 the following repre- 
sentations of the two algebras 

M 

K ±l ^i ±M J\a%, S^^r^c*, F^J2 lXc - K ^ ( 83 ) 

X—l X X 

and 

e x c x c* x+1 - clc x+1 +i(n x - n x+ i), x = 1,2, .... M — 1 . (84) 

One easily verifies that the above algebraic relation of U q (sl2) and TZ/m(<?) are satisfied. In 
addition we have the identities 

E 2 = F 2 = and K 2 = (-1) M . (85) 

Moreover, we have that the action of the Temperley-Lieb and the quantum algebra commute 
(this is a special case of the quantum analogue of Schur-Weyl duality [25]). 

[U q (sl 2 ),TL M (q)} = 0, (86) 

whence the Hamiltonian 

M-l 

.9 = 1: H = H'=J2 e * ( 87 ) 

x=l 

is quantum group invariant. For g = 1 these relations hold for both, odd and even numbers 
of sites. However, note that we have 

As we shall see below the vanishing of the commutator for M even is closely related to the 
fact that the discrete wave function of an associated fcrmionic zero mode has vanishing norm, 



This has profound consequences for the quasi-Hermiticity of the Hamiltonian, since the quasi- 
Hermiticity operator rj ceases to be positive definite, 77 > 0, and instead becomes positive 
semi-definite, rj > 0. One of the novel results in this article is the observation that a similar 
scenario also happens for M odd, albeit at the value g = ^-i > 1- ^ n order to discuss 
these issues and to interpret the new inner product in terms of representation theory it is 
favourable to introduce yet another algebra. 

Definition 3 The universal enveloping algebra [/(g[(l|l)) of the superalgebra gl(l|l) is the 

Ij-i-graded associative algebra (over C) generated by the odd (fermionic) elements {^ ± } and 
the even (bosonic) elements {Y, Z} satisfying the relations 

[Y,X ± ]=±X ± , [Z,Y] = [Z,X ± ]=0, [X+,X-]+=Z. (89) 

There is a natural anti-involution or ^-structure on U{q{{\\1)) by setting 

{X ± )*=X T , Y* =Y, Z* = Z . (90) 
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The following f/(0l(l|l) ^representation at g — 1 is a special case of the representation 
discussed in [TT] for the deformed case C/ ? (g[(l|l)) on which we will comment later, 



X+^Z^c:, X-~E^c„ Y^S*, ^{ M \ £™ • (91) 

Again, we observe a crucial difference for the two cases M even and odd: the vanishing of the 
central element Z. We will now discuss for g ^ 1 an extension of the above representation. 
For M even we shall see that (|91l) is a singular limit of a representation with Z^O. However, 
we start our discussion with the case M odd, since it is more well-behaved and shall then 
move on to M even pointing out similarities and differences for < g < 1. 

4.2 Odd number of sites 

For M odd the quantum group symmetry can be extended. As noted earlier we have for 
9 = ir/2 and M odd the following two Bethe roots z\ — ±V~1 giving rise to a zero mode of 
the Hamiltonian (we shall focus on z = i) 

[ff,c* =i ] = [#,4=i]=0, (92) 
where the corresponding wave function is given by 



™ — igcos ^ , 
^ z=l {x) = 2 2 . . (93) 

' M+l _ M-l 2 

2 2 y 



This solution prompts the following operator definitions 



U = J2sm^c x = £(-l)^c, (94) 



2 

x odd 



and 



such that 



and d z=i — — = ^ . (96) 



M+l M-l 2 / M+l M-l 2 
— 2 2~ ~~ ~ 9 

The operators t7, V satisfy the relations 

U 2 = V 2 = [U,V]+ = [U,V*]+ = [U*,V]+ = (97) 

and 

M + 1 M - 1 

CT]+ = — s~ ^ ]+ = — s~ 1 • (98) 



Thus, up to a trivial renormalization, U — » fZ/y (M + l)/2 and F — + V/ \J (M — l)/2, we can 
think of {/, V as two fermionic oscillators. Alternatively, we can also interpret them as two 
representations of the non-deformed superalgebra gl(l|l) introduced above by identifying 

Y^S Z , Zt^—£ll, X+^U*, X-^U (99) 

and 

M — 1 

Y^S Z , Z i — * — - — 1, X+^V*, X~^V. (100) 

We thus obtain two distinct representations of g[(l|l) which in addition preserve the *- 
structure (l9"01) with respect to the original inner product. For Q < g < 1 neither of these 
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two representations by itself give rise to a symmetry of the Hamiltonian, instead we have to 
consider the combined representation 

Y^S Z , Zh (M±l - g 2M-±^ i ; X + ^U*-igV*, X~ ^ U - igV . (101) 

According to ([5!?]) the generators X ± create respectively annihilate zero modes of the Hamil- 
tonian and we therefore have an £/(gI(l|l))-symmetry of the Hamiltonian (O, 



[# ff ,C/( t(l|l))]=O, Q< g< ^^±l, (102) 

Note that in the representation (|101[) the ^-structure (|90[) is not preserved with respect to 
the original inner product, 

X ± ? (X*)*, (103) 
but with respect to the new inner product induced by rj, 

r/X* = {X T )*ri . (104) 

This allows one to interpret the new Hilbert space structure induced by rj in a purely rep- 
resentation theoretic setting. It also singles out the J7(fll(l|l))-symmetry over the U q (sl2)- 
symmetry whose extension we discuss next. 

M+l 



Namely, for M e 2N + 1, < g < ^Jjfzj, and q = V^T we now set 

A-wfU b~ "--^ , f» zmm . (1 o5) 

x=1 \J— — 9 \J— — 9 

Employing the anticommutation relations for U, V and their Hermitian adjoints, one easily 
verifies that this representation is well-defined and that the relations (183)) continue to hold. 
Apparently, we recover the familiar representation (|83|) on the unit circle in the limit g — ► 1. 
The obvious guess for an extension of the Temperley-Lieb algebra for g > 0, 

e x c x c* +1 - c* x c x+ x + ig(n x - n x+1 ), x = 1,2, M - 1 , (106) 

does yield the modified commutation relations 

e l = (1 - 3 2 )[(! ~ n x )n x+1 + n x (l - n x+ i)] (107) 

and 

e x e x ±ie x = g 2 e x + ig(l - g 2 )(n x - n x+1 )[l + {n x±1 - n x+1± i)(n x - n x+1 )] . (108) 

This extension of the Temperley-Lieb algebra in terms of "local Hamiltonians" , i.e. the 
nearest neighbour terms, does not preserve the algebraic structure at g = 1. Furthermore, 
Schur-Weyl duality is broken for g < 1: the commutation relations between the fermionic 
oscillators and and the extended Temperley-Lieb generators (|106[) are 

[U,e x ] = { H^+i^c,), xodd 
[ (—) 2 (c x +igc x +i), x even 

iv,e x ] = { ;-£i c * +i5Cx+i j' " odd , (no) 

\ {-) 2 [c x +i - tgc x ), x even 
whence we now have for the quantum group generators the identities 



<>l-.</: JjfZJ : [E,e x + e x+1 ] = [F,e x + e x+1 ] = 0, x = 1,3,5, ...,M- 2 . (Ill) 
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Thus, while for M € 2N + 1 the Hamiltonian remains quantum group invariant for all values 



< .g < y' M-i ' we now have to consider pairs of the extended Temper ley-Lieb generators. 

Note that the representation (|105[) becomes singular at g 2 = (M + 1)/ (M — 1) , this is also 
the value where the fermionic modes X + = U* — igV* and X~ = U — igV anticommute, 



M + l 



[x + ,x-]+ = z = o 



(112) 



M —1 

This is precisely the scenario mentioned above for M even at g = 1. We see that the 
associated wavefunction (|93p becomes singular as its norm vanishes. Moreover, the Jordan 
normal form J of the Hamiltonian possesses now non-trivial 3x3 blocks. For instance, we 
have for M = 5 and the sector S z = 1/2 the Jordan normal form, 

/ 




J = 





1 







1 

: 











Some of the eigenvalues smoothly join up as can be seen in Figure 4. 

Spectrum of H g for M = 5 and S z = 1/2 




Figure 4. Spectrum of the Hamiltonian (ffj). 

Beyond the threshold value g 2 = |§zt the Hamiltonian has complex eigenvalues. Thus, 
quasi-Hermiticty of the Hamiltonian can only hold for g 2 < and this value approaches 
the unit circle as M — > oo. 
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4.3 Even number of sites 

For M even we define the operators U, V analogously to (|M|) and (|95p . and obtain again 
representations of g[(l|l). But we now have 

P, U*} + = [V, V*}+ = yl and Z ^^-(l-g 2 )l. (113) 

The definition of the other operators in pOip is unchanged. Due to the absence of a zero 
mode, we do not obtain as before a [7(fl[(l|l))-symmetry of the full Hamiltonian. Instead the 
generators X only commute with the truncated Hamiltonian where the last Temperley-Lieb 
generator is omitted, 

[Hf' nc ,C/( I(l|l))]=O (114) 

with 

M-2 A/-2 

Hf™ =H g - e M -i = £ e * = 2 E + « + i + <+i)] ■ (H5) 

x—1 x—1 

This already signals that this representation mimics the one of the chain with M — 1 sites. 
There is more evidence for this interpretation: for an extension of the £/ g (s^2 ^representation 
one now needs to introduce the operators 

M62N: K^=i^f\al E = ^=M=, p = _ ( ^+£j£ ; (116) 



which again satisfy the right commutation relations when setting g = \/— 1. Note, however, 
that the limit g — > 1 is now ill-defined and that we have replaced M — > M — 1 in the Cartan 
generator K. As for M odd we have the relations 

0<.g<l: [E, e x + e x+ i] = [F, e x + e x+ i] = 0, x = 1, 3, 5, M - 3 . (117) 

An alternative way to restore quantum group or J7(g[(l|l))-symmetry is by adding for 
M/2 even the generator 

ejw = cmc* - c^ci + ig(n M - "x) (H8) 
and consider the Hamiltonian with periodic boundary conditions 

M M 

Aperiodic = ^2 &x = y^( C a C x+X - 4 c x+l), M + 1 = 1 . (119) 

a:— X x—1 

The latter apparently does not depend on g and is simply the Hermitian Hamiltonian which 
describes free fermions on a lattice. We will not discuss this Hamiltonian further in the 
present context. 

Finally, for M even there is an alternative representation of the quantum group symmetry 
which is obtained from PT-reversal. Noting the transformation identities 

PUP={ ^"f^ M ° dd and PVP={ «T y ' M ° dd (120) 

| {-)~V, M even 1 {-)~U, M even v ' 

one easily verifies that 

^w^rU e i-i , - y-^ (i2i) 
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yields another representation of U q (sl2) with q = — \J — 1. The commutation relations ()1 17p 
are then modified accordingly by conjugation with the PT operator. 

Let us return to the case g = 1 for an even number of sites. If we only consider the 
fermionic oscillator modes X as before for M odd, then we see that they anticommute at 
g = 1, i.e. we have again Z = 0. We encounter the same scenario as before for M odd, the 
corresponding discrete wave function has zero norm, 

M 

M e 2N : ^CiW=0, 

x=l 

and the Hamiltonian has non-trivial Jordan blocks, however these are now of size 2x2. To 
be concrete we state here the Jordan normal form J of the Hamiltonian for M = 4 when 
restricted to the S z = subspace, 
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Looking at the spectrum of the Hamiltonian we see once more that some eigenvalues smoothly 
join up at g = 1, compare with Figure 5. 



Spectrum of H g for M = 4 and S z = 




0.0 0.2 0.4 0.6 0.8 1.0 

s 



Figure 5. Spectrum of the Hamiltonian (0. 

This signals the appearance of complex eigenvalues for \g\ > 1. Thus, quasi-Hermiticity only 
holds true inside the unit disc for M even. 

5 Perturbation theory on the imaginary axis 

While we have shown the existence of the map rj for boundary fields a — (3 inside the unit 
disc, its expression in terms of projectors on quasi-fcrmion states does not allow to easily 
compute the corresponding Hermitian Hamiltonian 

h = v^Hrj- 1 / 2 . 
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We therefore specialize now to the case < g < 1 and 9 = ir/2 in (|4"3")) and compute 
it perturbatively. We will arrive at the result already presented in the introduction: the 
presence of non-Hermitian boundary terms leads to a long range interaction in the bulk. 

Since r\ > we can write it in exponential form r\ = exp A. The quasi-Hermiticity relation 
then translates into 

H* = e A He~ A . (122) 

Employing the Baker-Campbell-Hausdorff formula the above identity can be expanded into 
multiple commutators 

H + y^ad h A H = H*, ad x Y=[X,Y], rj = e A . (123) 

fc>0 

This allows to determine the operator A perturbatively by separating the non-Hermitian 
Hamiltonian H into an Hcrmitian and non-Hermitian part of the following form, 

H = H Q + igH 1 , H Q = H£, Hi = H*, < g « 1 . (124) 

The Hamiltonian ([7]) has precisely this structure. Setting the boundary parameters in (fTJ) to 
purely imaginary values, a — —[3 = ig, g G R + , we define 

M-i z z 

H ° = -J2 t< c -+! ~ c - c -+i] ' Hl = aX ~° M - »i - n M (125) 

x=\ 

Since we have H* = H_ g we demand r\_ g = rj^ 1 leading to the following ansatz for a power 
series expansion of the metric operator in the coupling constant g, 

A = Y, 9 2n+1 Mn+i ■ (126) 

n>0 

From our previous considerations related to the appearance of Jordan blocks when the central 
element Z £ jj[(1|1) vanishes, we infer that the region of convergence for the above series is 
< g < 1 for M even and < g < y/(M + 1) / (M — 1) for M odd. 

Theorem 4 Collecting terms order by order in g we obtain the following equations for the 
coefficients in the expansion of the operator A, 

[H ,A 1 ]=2iH 1 (127) 

and for n > 1 

n 

[H Q ,A 2n+l ]^iY J ^k [A P1 ,...,[A P ^H 1 ]...} . (128) 

fe=l Pl + ---+P2fc=2n 

Here the sum runs over all compositions p = (pi, ■■■,P2k) of 2n with A Pi = if Pi is even. 
The coefficients are determined recursively via 

fc— 1 

Ax = 1/6 and X k = - g , * > 1 . (129) 

The first elements in this sequence are 

111 1 1 691 1 



v 6' 360' 15120' 604800' 23950080' 653837184000' 37362124800' 
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Proof. The recursion formula for the coefficients is proved via induction as follows. The 
case n = 1 is easily verified. Now assume that (If 28|1 holds true for all orders k — 1,2, ...,n. 
From the Baker-Campbell-Hausdorff formula 

k>0 

we find that 

n+l ^ 

[H Q ,A 2n+3 ] = Yl [ A Pi,~dA>^ H l\~] 

fc=l ^ pi + ...+p 2 fc=2n+2 

n+l 

+ X1 C2jfc + lV 51 [ A Pl,-,[ A P 2k+l, H o]-] ■ 

k=l K >' pi+...+p 2 fc + i=2n+3 

Employing our assumption we can rewrite the second term on the right hand side in terms 
of Hi. For P2k+i = 1 we find 

[A>2s>+i > Ho] = -2iH\ 
and for f < p2k+i < 2[n + 1 — k) + 1 we have 

P2fc + 1~ 1 
2 

[A^ +1 ,#o] = -« X i E [A qi ,...,[A q2j ,Hx]...] . 

i=l gi+...+<j2j=P2fc+i-l 
Inserting these expressions, rearranging sums and collecting terms we find 

71+1 f 1 9 k- 1 A 1 

[H 0! A 2n+3] =iJ2\ m - ^nyy - £ (^Syy E , ^ , *].••] 

fell 1 ' V 7 3-1 W 1 7 J Pl +...+ P 2fc=2n+2 

which yields the desired recursion formula for the coefficients. I 

In order to facilitate the comparison with results in the literature, we explicitly state the 
identities up to order seven: 

[H , A x ] = 2%H X 

[H ,A 3 ] = Ua^A^Hx]] 
6 

[H ,A 5 ] = ^[Ax,[A 3 ,Hx]] + ~[A 3 ,[Ax,Hx]]-^[Ax,[AxMiMi,Hi]]]] (130) 
6 6 360 



and 



[H ,A 7 ] = l -[A 3 ,[A 3 ,Hx]] + ^[Ax,[A 5 ,Hx]]+ l -[A 5 ,[Ax,Hx] 
6 6 6 

Pl+P2+P3+P4=6 

[AiJAxJAxJAiJAi,^,^]]]]]] 



15120 



Note that our general formula differs from equation (3.20) in [T7] but reproduces the identities 
in equation (34) of [It)] . 

One can now use the identities (| 128[) to compute the metric operator A order by order. 
In addition to the latter equations one also easily deduces from the properties rj ought to 
obey that we must have 

[A, S z ] =0, PA = -AP and A* = A . (131) 
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Define 

at, y = c* x c y ± c x c* y , x < y (132) 

then we find that up to order 1 1 the terms in the expansion of the matrix A are of the general 
form 

. M— 2n— 1 ri— In— p 



A2n+1 — rr ^ a t,x+2n+l + * K i™' P ^ ( a x,a;+2p+l + a ll-x-2p,M+l-x) • (133) 

^ ' x = l p=0 x=l 

Explicitly they read, 

M-l 

A i = i Yl a ^+i 

X=l 

1 1 

^3 = g- ^2 a Xx+3 ~ qI ( a 1^2 + a M-l,Mj 
x=l 

■ M-5 . . j - . 

^5 = g a x,x+5 + 7^ ( a t2 + a M-l,Af) ^Y20 ( fl ^ 3 °M-2,M-l) _ J20 a ^-f-3,Af) 



x=l 



.4 



M-7 

,+ 

e+7 



+ 

M-3.M-2 



= 7i Yl a t* 

x — l 

~ 24to (°^ 2 + a ^-L^) ~ 48* + a "-W-i) ~ 840i (°^ 4 + ° 

+ 60i (^ 4 + °^-3- M ) ~ 560z + a ^-4,M-i) - (<6 + a M- 5 ,Ai) 

and 

. M-9 

^9 = g Y^ a x,x+9 + M ( a l,2 + a M-l,Af) + ( a 2,3 + a A/-2,Af-l) 
x — l 

+ ( fl 3,4 + a M-3,M-2^) + J^M ( a 4,5 + a M-4,M-3J + + a M-3,M J 

+ ilH>I ( a 2 : 5 + a A/-4,A/-l) + WPIOi ( a 3,6 + a M-5,M-2 



13440 ( fl l,6 + a A 



59 'ao, + a7 ^ 1 1823 



~ 13440 V "1,6 ~ "M-5.M ^ ^ 8064i ^"2,7 ~ "M-6,M-1 J ~ 40320i \"1,8 ~ "M-7,M 

The stated solutions can be checked by employing the relations 
Wo, a x.v] = a L,+i - .v + af 



J'"x,i/J ~~ u x,y+l "x— 1,3/ ~ "x,j|-l "x+l,y 

[-ffi.^x,!/] = (<W + s y,M)a* y , x<y 

and 

[ a a;,j/i a r,s] = ^y,r a x,s ^y,s a x,r ~ ^x,sO- r ,y ~ ^x,rd s ,y ■ 

Note that some of the terms in the general formula (|133p vanish if 2n + 1 > M. The general 
solution still needs to be found. Once the metric operator A is obtained, the corresponding 
Hermitian Hamiltonian h can be computed order by order. Define the Hcrmitian Hamiltonian 
h according to 

oo 

h = e*He-$ =H a + Y / 9 2n h2n , (134) 

71=1 

where the same region of convergence is implied as in the case of A. 
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Theorem 5 The terms in the series expansion of h are given by 

n 

h 2n = iJ2K E [Ap,,...,^,^]...}, (135) 

fc=l pi+...+p 2 fc-l=2n-l 

where the sum runs over all compositions p — (jc>i, ...,j>2fc-i) of 2n — 1 and the coefficients 
\' k are computed from the coefficients \k in 11128]) via the formula 

, = 2k - 1 _ ^ Xk-j n „ fi x 

k 2 2fc ~ 1 (2fc)! 2^{2j)\ ' 1 ' 

The first terms in the above sequence are 

,1 1 1 17 31 691 

~ ( 4'~i~92' 7680 '~~ 5160960' 371589120' 326998425600' 

Proof. An induction proof similar to the previous one for the series expansion of the 
operator A. ■ 

Again, we state the first terms explicitly in order to allow for comparison with the liter- 
ature (see equation (3.24) in [17]), 



h 2 = 


i 

4 


[Ai,Hx] 




hi = 


% 
4 




i 

~ 192 


h 6 


i 
4 


[A^H{\ 


i 

+ 192 



E [^,[^,[^3,^1] 

Pl+P2+P3=5 

' [Aufa^AulAiAMiHi]]]]]- 



7680 

Note that h only depends on g 2 as it must due to H* = H- g . Inserting the explicit expressions 
for the series expansion of A up to order n = 5 we find the Hamiltonian stated in the 
introduction. 



6 Exact results for short spin-chains 

In this section we present exact expressions for the quasi-Hermiticity operator rj and the 
operator A employing the fact that for small M there are only a few terms appearing in the 
expressions for the operator A, compare with (|133|) . That is, motivated by our perturbation 
theory results we make the ansatz 

[M/2-1J M-2p-l 

^ = E E ^x^ ( a t,x+2p+l + a M-x-2p,AI+l-x) ■ (137) 

p— x— 1 

For M small enough one can exponentiate the resulting matrix and solve the intertwining 
relation in order to determine the coefficients Xx . From our perturbative computation wc 
can also make an educated guess about the form of the Hermitian Hamiltonian h in order to 
facilitate the computation of the similarity transformation. Namely, we make the ansatz 

[Af/2-lJ M-2p-l 

^ = E E £x P \ a x,x+2p+l + a M-x-2p,M+l-x) ■ (138) 

p— x—1 
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Here we have used that h is parity invariant. This follows from the identities 

Prfi = p e ^/ 2 = e ~ A/2 P = jf^P 

and 

Ph = Prj^Hr}-? =r)-iH*r]ip = hP . (139) 

In the last part of this section we present another exact approach which is valid only 
at g = 1 for M odd but it exploits the graphical calculus attached to the Temperley-Lieb 
algebra and is interesting in its own right. 



6.1 Example M = 3 

It might be instructive to start out with the easiest, albeit not the most interesting, case 
M = 3. The non-Hermitian Hamiltonian reads 



H 



n,2 



a 2,3 + *5( n l - n 3) 



(140) 



where a x y has been defined in (|132j) . The operator A ~ In 77 is then computed to be 



.4 



24- q 

arccos 

^ 9 



V2 



(at 



l 2,3) 



(141) 



After exponentiation we obtain the quasi-Hermiticity operator and from its explicit form, 



V = 



( 1 







V 







9 ~ 
2ig 

g 2 -2 

V 










2ig 
9 2 -2 
-1 - ■ 
2»ff 







2 -2 
2ig 

2 

-g^2 



9 2 -2 
2% 



g 2 -2 
r 



9" 





2»g 

-1 - - 

9 

2ig 



9-1 
_ 2ig 







S 2 -2 
2 

"F^2 



\ 







1 / 



(142) 



one can explicitly verify that 77 satisfies all the necessary conditions. Moreover, as predicted 
from our algebraic analysis related to the J7(g[(l|l))-symmetry of the Hamiltonian, we see 
that 77 has a pole precisely at the value g 2 = (M + l)/(Af — 1) = 2 where the central element 
Z e J7(gt(l|l)) vanishes. The renormalized 77, 

77 -> (g 2 - 2)77', 

still intertwines the Hamiltonian, rj'H = H*r( at g 2 = 2, but is not any longer positive 
definite but semi-definite. In principle one can also compute the square root for 77 1 / 2 for 
general values of g but as the expressions become rather unwieldy we restrict ourselves to 
stating the result for g = 1, 



lim 77 2 = 
9^1 



/ 1 


























\ 





1 

2 


+ 71 


i 


1 1 

2 ^2 























i 




i 





















1 

2 




i 

V2 


1 + 1 

2 ^ V2 
































1 

2 


+ 7f 


i 

V2 


1 

2 


1 

V2 





















i 

75 1 




i 





















1 

2 




i 


2 > 


1 





V 


























1 / 



(143) 
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After performing the similarity transformation we end up with the Hcrmitian Hamiltonian 



lim h — 



a l,2 + a 2.3 

V2 



(144) 



Arguably, we do not see here the manifestation of the long-range nature of the interaction, 
since the spin-chain is simply too short. However, this case illustrates the working of the 
general ideas and concepts introduced earlier. 



6.2 Example M = 4 

We saw previously that the cases M even and odd are different and it is therefore important 
to have another simple but concrete example. We now hnd 



-( a i2 + a U) + 



1 2,3 



+ 



1,4 



(145) 



where the coefficients are given by 
£ = 2 arctanh ■ 



1 



r 



( = 2i arctanh ■ 



-1 + VT 



_ (146) 

Inserting these expressions into the formula for A does not lead to any further drastic sim- 
plifications. Also the expression for the quasi-Hermiticity operator in the spin-sector S z = 
looks now more complicated: after rescaling 



7 



rf = 7(3 - g 2 - 2 7 )?7 



one finds the matrix 



1 


ig 


1-9 2 - 


7 


l-3 2 -7 


ig (2 - g 2 - 2 7 ) 


5 2 + 2 7 


- 2 


* 


3 - 2 7 


-ig (7 - 
7 - g 2 (7 


2) 


-*S (7 -2) 


~g 2 


ig (2 - ff 2 


-27) 


* 


* 


-2) 


2-27 


-iff (7 -2) 


1-5 2 - 


- 7 


* 


* 


* 




7 - g 2 (7 - 2) 


-iff (7 -2) 


l- ff 2 - 


"7 


* 


* 


* 




* 


3-2 7 


iff 




* 


* 


* 




* 


* 


1 





y** * * * 1 j 

(147) 

Here we have omitted the lower diagonal part, since rf — 77. Again one can explicitly check 
that this operator satisfies all the necessary requirements. From our previous discussion we 
expect that in the limit g — > 1 the operator 7/ should cease to be positive definite and become 
semi-definite. In fact, one verifies that this is the case by using the expression 



lim rj 

9-1 



1 


i 








i 


-1 


— i 


3 


2i 


2i 


-1 


i 





-2i 


2 


2 


2i 








-2i 


2 


2 


2i 





— i 


-1 


-2i 


-2i 


3 


i 


-1 


—i 








— i 


1 



(148) 



6.3 Example M = 5 

Finally to present one example where the long range nature at g = 1 is slightly more apparent 
we consider the M = 5 chain. We specialize from the start to the case g = 1 and find that 



lim A 



(«1,2 



^4.5/ 



72 



-(a 



2,3 



2 3.4) + : o (°1,4 



*2,5 



(149) 



27 



with 



i In 5 



C = i In 



2 + ^ + 2 v« 



'4 + V5 



After exponentiation we obtain the quasi- Hermiticity operator. In the spin S z = 3/2 sector 
it reads explicitly 



lim T7 = 



iV5 



2 + 



V 



V 5 

1 + V5 

5+3V5 
5i 

-1 - 



-1 3_ _J_ _j 

I 3* 1 2 



1 



v/5 



1 7= — ^= — 2 

V5 iV5 



4 



3i 



V5 
o±3V§ 
5i 



tT V5 

1 + V5 



V5 



-1 L. _a o 

V5 iV5 



-1 

i + 
1 + 



and its square root is computed to be 



3i 
V5 



(150) 



lim 77 2 

9->l 



7+9 ^5-^+-!^. 


ity5 


l_ 5 ^5+3V5-9/5 1 / 4 


— i 


7-2 4 / 5- v / 5+-A 4 - 

2 . 


22 




22 


2^5 


22 


2 

1-5^5+3^-^ 


2 

— i 


i 

4+2y5+v / 5+8/5 1 / 4 


— 1 

2^5 


i 

1-5^5+3^-^4 


22 


^5 


11 


^5 


22 


i 

2^5 

7-2^5-^+-^ 


-1 

2^5 


— i 


2 

^5 


2 

7+9^5-^+^4 


22 


2^5 


22 


2i 


22 



av ; ; v ■jt^v -j— g/ yo 5 1 /" 

\ 22 2^5 22 2i 22 / 

(151) 

Using these matrix expressions one can compute the Hermitian Hamiltonian h and we find 
that it is of the expected form, 



lim h = p x (a X)2 + a 4 5 ) + p 2 (a 2 3 + a 3 4 ) + p 3 (a 1)4 + a 2 5 ) 



(152) 



with 



_ 9-6^-^/2(15+23^) _ 3-275-V40+21y^ _ -2+5^-^2(15+23^5) 

Pi — to ) ^2 — 11 > P3 — 22 ' y 106 ) 



The non-local nature of the interaction is now visible in terms of the nonvanishing coefficient 
p 3 . Below we show the result for the spin sector S z = 1/2, however we compute it via 
different means. 



6.4 The quasi-Hermiticity operator and the Temperley-Lieb algebra 

In this section we indicate an alternative way to compute the quasi-Hermiticity operator rj at 
9=1 when M is odd employing the Temperley-Lieb algebra. From the expressions (jT4")l we 
expect that in general rj will be highly non-local in terms of the spin basis, since the quasi- 
particle creation and annihilation operators involve sums over the argument of the discrete 
wave functions. Thus, if we compute the matrix elements of rj in a fixed spin sector then all 
elements will be non- vanishing. Our example above for M = 5 and S z = 3/2 confirms this, 
compare with (| 1 50[) . However, there is another choice of basis in which rj simplifies, that is 
where many of its matrix elements are zero. Namely, we consider the q — > i = \f— 1 limit 
of the dual canonical basis [26, 27J which we denote by {U}. The elements of this basis are 
in one-to-one correspondence with algebra elements a% € TLm and the latter can be written 
down in terms of Young tableaux, see e.g. [2"5] . 
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Fix a spin sector, S z — const., and set m = M/2 — S z . Let A„ be the rectangular Young 
diagram with n rows of TV — n boxes, 



> to 



M—m 



Then we assign to each subdiagram A' C A m a vector as follows. Let t be the unique standard 
tableau (column and row strict) of shape A' whose entries are consecutive integers with entry 
n in the upper left corner. For example, 



m 


m + 1 


m + 2 




s 


m — 1 


m 




s - 2 














s> 











, to < s < M, 1 < s' < to . (154) 



Reading the entries of the tableau from left to right and top to bottom we set 

1 1 > e s /e s '_i • • • f-s~2 ' ' ' 6 m — ie s • • • e TO _|_ie m f2 m , (155) 

where 

TO = u_ <g> v_ • • • ® w_ (8) v + <g> u + • • • <g> w+ (156) 

S v ' 

m 

is the vector corresponding to A' = 0. Note that for fixed to there are as many of these 

tableaux as the dimension of the spin sector, namely ( ) . 

\m J 

Example. Let M = 5 and m = 2 then we have the following Young diagrams and tableaux: 
t = 0, [2l, 



2 3 



2 


3| 







2 


3 


1 


2 



2 3 4 



2 


3 | 4 | 


2 


3 


4 


2 


3 


4 


1 


7 


1 


2 


7 


1 


2 


3 



The corresponding algebra elements a £ TLm(<z) are 

a = 1, e 2 , eie 2 , e 3 e2, e4e3e2, eie3e2, e 2 eie3e 2 , eie4e 3 e 2 , e 2 eie4e3e2, e3e 2 eie4e3e2 . 

Each of these algebra elements we can also represent as a link or Kauffman diagram. Iden- 
tifying 

<?,- = 



i i+1 

and realizing multiplication by concatenation from above, we find the following diagrams for 
the algebra elements, 



a £ TL b 


Oi = 1 


02 = e 2 


«3 = eie 2 


a-4 = e 3 e 2 


a 5 = e 4 e 3 e 2 


«6 = eie 3 e 2 


































diagram 






/ 











and 
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a 6 TL 5 


aj = e 2 eie 3 e 2 


a s = eie 4 e 3 e 2 


09 = e 2 eie 4 e 3 e 2 


aio = e 3 e2eie4e 3 e2 








yy \ yy 






diagram 













We will make use of these diagrams momentarily. Employing the representation TLm — * 
End V® M introduced earlier, 

e x !-» c x c* x+1 - c*c x+ i + i(n x - n x+ i), 

we can now generate the basis vectors U by acting with each corresponding algebra element 
di onto the vector f2 m . Computing from our previous expression (|T4|) the quasi-Hermiticity 
operator rj we can evaluate its matrix elements in this new basis. Based on numerical com- 
putations for M = 3, 5, 7 we arrive at the following conjecture. 

Conjecture 6 Denote by G the Gram matrix of the dual canonical basis vectors {ti} with 
respect to the rj-product, i.e. 

Gij — {tii V^j) ■ 

Then we have 

Gij = whenever tr(ai<2j) — 0mod2, 
where a,i,aj are the algebra elements corresponding to U,tj and 

tr a = tr ocm = number of closed loops 

which are obtained by closing the planar diagram associated with a. An example for a = e 2 
is shown below, 




= M- 1 . 



Remark. Note that the above relation does not necessarily imply that Gy ^ if tr(aiOj) = 
1 mod 2. The Gram matrix inherits from r\ the further properties 

det G = 1, G > and G M = G 3i £ K . 

The knowledge of the Gram matrix is sufficient to compute matrix elements with respect 
to the 77-product. In particular the intertwining property, r\H = H*rj, translates into the 
following identity for the Gram matrix 

GH = U l G (157) 
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where the matrix Ti. has only integer entries, Tlji G Z, and is defined via 



Hti - ^ 1 : 



(158) 



Note that Ti. can be computed graphically using link diagrams. For instance, we find for our 
example M — 5 and S z = 1/2 stated above that 
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The corresponding Gram matrix is 



G 



f 2(3+x/5) 
5 



2(l+V5) 
5 

1 , J_ 

5 ^ V5 





3 



1 



2(1+75) 
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2(2+^5) 

5 
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3 
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5 
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75 
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75 







3(3+ Vs) 
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I 

5 



9 
5 





1 + 





1 





2 



f 



1 

f 



f 



(159) 



(160) 



One can check that the above conjecture holds true for this example using the stated link 
diagrams. 

In addition to the identity (|157[) we obtain further constraints on the Gram matrix by 
employing PT-invariance, 77" 1 = PrjP = fj = rf, from which we conclude that 



G, 



and hence 



M*GM = G, PTU = h M 3 



(161) 



It is desirable to find a closed formula for the Gram matrix elements similar to the case 
treated in fl9l [201 . 



7 Results for boundary fields off the imaginary axis 

In this section we briefly discuss some aspects when a = (3 do not lie on the imaginary axis. 
This case for a — (3 on the unit circle has been investigated previously pTJ [311 [12] albeit not 
in the context of quasi-Hermiticity. Here we relate our discussion to these previous results. 
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7.1 UJgl(l\l)) invariance and Hecke algebra 



In order to discuss the quantum group symmetry we consider the Hamiltonian H' instead of 
H, 



A I 



H' — — 



T m°~m+l + <7 m°'rn+l 



'm+l 



m—l 



This Hamiltonian can be viewed as an element in the Hecke Algebra, compare with the 
discussion in 



Definition 4 The Hecke algebra H. n (q) with q G C is the associative algebra (over C) ob- 
tained from the generators b\, ...,b n —i subject to the relations 



b i b i = b t h = 1, bjbj+ibj = bj + ibjbj + i, b z bj = bjbi for \i - j\ > 1 



and 



b 2 i +{q-q- 1 )b i = l. 



Setting n — M and q = —a 1 one easily verifies that 
bi i ► CiC* +1 - c*c i+1 - cT x ni - a(n i+1 - I) 



/ q \ 

q-q- 1 1 



and 



b i 1 h-> CiC* +1 - c*c i+ i - a 1 (n i - I) - an i+1 



10 
V -q-i J 



( -q- 1 \ 

10 

I q- x -q 

V 



yield a representation of the Hecke algebra 7Ym(<Z = x ) over V® and, furthermore, we 
have the identity 



b t + b l 



(162) 



There is a action of the quantum group U q (gl(l\l)) [5S1[3D] which is "dual" to the action of 
the Hecke algebra in the sense that both of them commute, i.e. the Hamiltonian is quantum 
group invariant. 

Definition 5 Consider the associative algebra U q (Ql(l\l)) over C(q) generated by the ele- 
ments {X ± ,y ±1 , Z ±1 } subject to the relations 

zz- 1 = z~ x z = yy~ x = y~ x y = 1, yx ± y- 1 = q ±1 x ± , [y,z] = [z,x±\ = o (163) 

and 

'± v±l — n r v + r i — ~ ~ (164) 



[X ± ,X ± ]+ = 0, [X+,X-]+ = Z Z _. 

q-q 1 



Setting q = —a 1 and identifying 



Z ±X ^q ±M , 



, ,4- \ M + l 



\ ^ M+l 



we obtain a representation of this algebra over V® M . Moreover, its action commutes with 
the action of the Hecke algebra, 



[U g (gl(l\l)),H M (q)}=0 ) 
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and, thus, with the Hamiltonian H' . As before the quantum group invariance is reflected 
in the existence of a particularly simple solution for the discrete wave function. Setting as 
before q = —oT l = e 10 the mentioned wave function and creation operator are 

x — 1 

Using these expressions one finds 

[H,c* g ] = -2cos6> c* e and [H',c* g ] = . (166) 

Similar as in the case 9 = it/ 2 and M even the Hamiltonian has non-trivial Jordan blocks 
when the norm of the wave function becomes singular, that is if we choose 9 such that 

sin MB _ 
sin# 

Also the representation theoretic interpretation of r\ extends from 9 = iir/2 to the present, 
general case: the new inner product (•, 77-) gives rise to a representation of the quantum group 
£/g({j[(l|l)) with the correct ^-structure, 

7]X+ = (X-)*r) . (167) 

The above identity has another consequence. Since we also have that PX + P — (X~)* it 
follows that the C-operator obeys 

[C,Z4( [(1|1))] = O, C = P V 

and hence must be an element in the Heckc algebra. This is precisely the line of reasoning 
employed in [9] to find an algebraic expression for C. However, we need an additional 
ingredient to complete the analogous computation: the decomposition of the state space into 
Wq(0[(l|l))-modules [32]. We leave this problem to future work. 

7.2 Perturbation theory for arbitrary 8 

In this section we wish to confirm that our previous picture for 9 — n/2 stays intact for 
general 9. Namely, we wish to show that also here the non-Hermitian boundary fields in 
H, H' correspond to non-local hopping terms in h. 

For arbitrary 9 we have to modify our previous approach to the perturbative computation 
of 77 = exp A. We now define 

a = 13* = -ge l \ < g < 1 : H g = H + gB x 

with 

M~ 1 id z 1 —id 7 

H ° = - E [< c -+! ~ c -<+i] » H ^ = ~ - ° o ° M = -( e ^ ni + e " iS "M) • (168) 

x=l 1 

For < g <C 1 we again expand the operator 77 = exp A using the Baker-Campbell-Hausdorff 
formula but now we cannot make the ansatz that A only depends on odd powers of the 
coupling constant g, since H* ^ H— g . Thus all powers of g are occurring in the series 
expansion, 

A = Y J 9 n An- (169) 

n>0 
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Collecting once more terms of the same order in g we now arrive at 



[H , At] = 2H_ 

[H ,A 2 ] = [A 1 ,H 1 ] + ^[A 1 ,[A 1 ,H ]} = [A 1 ,H + ] 

[H , As] = [A 2 ,H + ]-±-[A l ,[A l ,[A l ,H Q ]]] = [A 2 ,H + ] + \[A l ,[A l ,H-]l (170) 

12 



where 



H + = - — - = —cos8(ni +um) and H- = — - = ism.6{nM — n\) . (171) 

The explicit results for the first two terms in the expansion of the matrix A are, 

M-l 



M-2 

A 2 = Yl [(l + e 2ie )c*c x+2 -(l + e- 2i Vc: +2 ] 



The corresponding Hermitian Hamiltonian h now reads up to third order in the coupling g, 



with 



hi = H + = — cos 6{n\ + nj^), 

h 2 = 7 [Ai,ff_] = ^-e^(ctc 2 + c^_ 1 c M )+h.c., 
4 4i 

^3 = J [A 2 ,H-] = ^-(l + e l2e )(c* 1 c 3 + c* M _ 2 c M )+h.c. . 

Here "h.c." stands for the Hermitian conjugate of the previous term. From this we infer 
that for general values of 9 also interactions between sites separated by an even number are 
possible. 

8 Conclusions 

In this article we encountered numerous new aspects of the XX spin-chain with non-Hermitian 
boundary fields. While this simple integrable model has been studied intensively before in 
the literature, it has recently received renewed attention in connection with logarithmic field 
theories where the main focus is on representations of the Temperley-Lieb algebra with non- 
trivial Jordan blocks [71 [S] . 

The present work has highlighted a very different aspect: purely on physical grounds one 
wishes to have a Hermitian quantum Hamiltonian (without Jordan blocks) in order to ensure 
a unitary time evolution of the system. This is one of the essential demands of quantum 
mechanics. We have seen that this can be achieved for complex boundary fields with values 
inside the unit disc via two different, albeit closely related routes: one can either introduce 
a new inner product or perform a similarity transformation to a Hermitian Hamiltonian. 
The latter leads to a new physical interpretation of the non-Hermitian XX spin-chain: it 
corresponds to a free fermion system with long range hopping. The probability of long range 
hopping taking place is controlled by the absolute value of the complex boundary fields. This 
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new perspective on the XX spin-chain with non-Hermitian boundary fields entails a range of 
other physically interesting questions, such as finite size effects and correlation functions. We 
already touched upon the finite size scaling of the groundstate energy in the text, since the 
latter contains information about the respective CFTs in the thermodynamic limit and thus 
would connect with the discussion in [6l [TJ [8] . Due to the absence of an exact solution for 
the Bethe roots when the boundary fields lie within the unit disc, we were unable to obtain 
conclusive results. An alternative approach might be to find a field theoretic model which 
allows one to compute the partition function, similar as it has been the case for critical dense 
polymers on the lattice [BJ. 

To find the correct field theoretic counterpart might be facilitated by the algebraic struc- 
tures pointed out in this article. We explicitly constructed representations of the quantum 
group Uq{si2) with q = and the superalgebra J7(fj[(l|l)) which can be extended from 

the unit circle along the imaginary axis. Any field theory describing the thermodynamic 
limit should reflect these algebraic features. Recall that for M odd these algebras provided 
symmetries, while they did not for M even. Our careful analysis showed that we could tic 
the appearance of (non-trivial) Jordan blocks in the Hamiltonian to the vanishing of the 
central element among the generators of U(gl(l\l)). For M even this happens precisely if 
the boundary fields lie on the unit circle, but the new result here is that it also happens for 
M odd just outside the unit circle. Moreover, we recall from the main text that for M even 
one obtains 2x2 Jordan blocks, while for M odd one has 3x3 blocks. Thus the algebraic 
picture put forward in (8] for M even needs to be modified for M odd and an interesting 
problem is whether one can find also here corresponding logarithmic CFTs. 

The superalgebra has also been connected to the discussion of quasi-Hermiticity. 

The new inner product preserves the natural *-involution of U(gl(l |1 )), thus providing a rep- 
resentation theoretic interpretation of our construction, ft is this feature which singles out 
non-Hermitian quantum integrable systems: due to their underlying algebraic structures they 
allow for various interpretations and exact constructions of the quasi-Hermiticity operator. 
This connection between abstract mathematical structures and physically motivated concepts 
such as PT-symmetry is mutually beneficial, it leads to new physical insight and interesting 
mathematical questions even for a simple and well-studied model such as the XX spin-chain. 
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